REPORT  DOCUMENTATION  PAGE 


Form  Approved  OMB  NO.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions, 
searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments 

regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including  suggesstions  for  reducing  this  burden,  to  Washington 

Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington  VA,  22202-4302. 
Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  oenalty  for  failing  to  comply  with  a  collection  of 
information  if  it  does  not  display  a  currently  valid  OMB  control  number. 


PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY) 

2.  REPORT  TYPE 

New  Reprint 

3.  DATES  COVERED  (From  -  To) 

4.  TITLE  AND  SUBTITLE 

Switching  Exponent  Scaling  near  Bifurcation  Points  for 

Non-Gaussian  Noise 

5a.  CONTRACT  NUMBER 

W91  INF-08-1-0336 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

411359 

6.  AUTHORS 

L.  Billings,  I.B.  Schwartz,  M.  McCrary,  A.N.  Korotkov,  M.I.  Dykman 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

- 1 - 

7.  PERFORMING  ORGANIZATION  NAMES  AND  ADDRE: 

University  of  California  -  Santa  Barbara 
Office  of  Research 

The  Regents  of  the  University  of  California 
Santa  Barbara,  CA  93106  -2050 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND 
ADDRESS(ES) 

U.S.  Army  Research  Office 
P.O.Box  12211 

Research  Triangle  Park,  NC  27709-2211 

12.  DISTRIBUTION  AVAILIBILITY  STATEMENT 
Approved  for  public  release;  federal  purpose  rights 

13.  SUPPLEMENTARY  NOTES 
The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and  should  not  contrued  as  an  official  Department 
of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other  documentation. 

14.  ABSTRACT 

We  study  noise-induced  switching  of  a  system  close  to  bifurcation  parameter  values  where  the  number  of  stable 
states  changes.  For  non-Gaussian  noise,  the  switching  exponent,  which  gives  the  logarithm  of  the  switching  rate, 
displays  a  non-power-law  dependence  on  the  distance  to  the  bifurcation  point.  This 

dependence  is  found  for  Poisson  noise.  Even  weak  additional  Gaussian  noise  dominates  switching  sufficiently  close 
to  the  bifurcation  point,  leading  to  a  crossover  in  the  behavior  of  the  switching  exponent. 

15.  SUBJECT  TERMS 
noise,  non-Gaussian,  bifurcation 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

ABSTRACT 

UU 

15.  NUMBER 

OF  PAGES 

19a.  NAME  OF  RESPONSIBLE  PERSON 

John  Martinis 

a.  REPORT 

UU 

b.  ABSTRACT 

UU 

c.  THIS  PAGE 

UU 

19b.  TELEPHONE  NUMBER 

805-893-3910 

8.  PERFORMING  ORGANIZATION  REPORT 
NUMBER 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 
ARO 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

54568-PH-OC.5 


Standard  Form  298  (Rev  8/98) 


Prescribed  by  ANSI  Std.  Z39. 1 8 


Report  Title 

Switching  Exponent  Scaling  near  Bifurcation  Points  for  Non-Gaussian  Noise 

ABSTRACT 

We  study  noise-induced  switching  of  a  system  close  to  bifurcation  parameter  values  where  the  number  of  stable  states 
changes.  For  non-Gaussian  noise,  the  switching  exponent,  which  gives  the  logarithm  of  the  switching  rate,  displays  a 
non-power-law  dependence  on  the  distance  to  the  bifurcation  point.  This 

dependence  is  found  for  Poisson  noise.  Even  weak  additional  Gaussian  noise  dominates  switching  sufficiently  close  to 
the  bifurcation  point,  leading  to  a  crossover  in  the  behavior  of  the  switching  exponent. 


REPORT  DOCUMENTATION  PAGE  (SF298) 
(Continuation  Sheet) 


Continuation  for  Block  13 


ARO  Report  Number  54568. 5-PH-QC 

Switching  Exponent  Scaling  near  Bifurcation  Po 


Block  13:  Supplementary  Note 

©2010  APS.  Published  in  Physical  Review  Letters,  Vol.  104,140601,  (2010),  (0601).  DoD  Components  reserve  a  royalty-free, 
nonexclusive  and  irrevocable  right  to  reproduce,  publish,  or  otherwise  use  the  work  for  Federal  purposes,  and  to  authroize 
others  to  do  so  (DODGARS  §32.36).  The  views,  opinions  and/or  findings  contained  in  this  report  are  those  of  the  author(s)  and 
should  not  be  construed  as  an  official  Department  of  the  Army  position,  policy  or  decision,  unless  so  designated  by  other 
documentation. 

Approved  for  public  release;  federal  purpose  rights 


PRL  104,  140601  (2010) 


PHYSICAL  REVIEW  LETTERS 


week  ending 
9  APRIL  2010 


Switching  Exponent  Scaling  near  Bifurcation  Points  for  Non-Gaussian  Noise 

Lora  Billings,1  Ira  B.  Schwartz,2  Marie  McCrary,1  A.  N.  Korotkov,3  and  M.I.  Dykman4’* 

1  Department  of  Mathematical  Sciences,  Montclair  State  University,  Mon  tclair,  New  Jersey  07043,  USA 
2U.S.  Naval  Research  Laboratory,  Code  6792,  Nonlinear  System  Dynamics  Section,  Plasma  Physics  Division, 

Washington,  D.C.  20375,  USA 

3 Department  of  Electrical  Engineering,  University  of  California,  Riverside,  California  92521,  USA 
4 Department  of  Physics  and  Astronomy,  Michigan  State  University,  East  Lansing,  Michigan  48824,  USA 
(Received  22  September  2009;  published  7  April  2010) 

We  study  noise-induced  switching  of  a  system  close  to  bifurcation  parameter  values  where  the  number 
of  stable  states  changes.  For  non-Gaussian  noise,  the  switching  exponent,  which  gives  the  logarithm  of  the 
switching  rate,  displays  a  non-power-law  dependence  on  the  distance  to  the  bifurcation  point.  This 
dependence  is  found  for  Poisson  noise.  Even  weak  additional  Gaussian  noise  dominates  switching 
sufficiently  close  to  the  bifurcation  point,  leading  to  a  crossover  in  the  behavior  of  the  switching  exponent. 

DOI:  10.1103/PhysRevLett.  104. 140601  PACS  numbers:  05.40.Ca,  05.70.Ln,  72.70.+m 


Physical  systems  display  generic  features  near  bifurca¬ 
tion  parameter  values  where  the  number  of  stable  states 
changes.  In  this  range  the  dynamics  is  controlled  by  a  slow 
variable,  a  soft  mode.  Its  noise-induced  fluctuations  are 
comparatively  large.  They  ultimately  lead  to  switching  of 
the  system  from  the  stable  state.  Close  to  a  bifurcation 
point  the  switching  rate  W  becomes  appreciable  even 
where  far  from  this  point  it  is  exceedingly  small,  for  a 
given  noise  level  (for  example,  for  given  temperature).  The 
high  sensitivity  of  the  rate  to  the  system  parameters  has 
been  broadly  used  to  determine  parameters  of  Josephson 
junctions  and  Josephson  junction  based  systems  [1-4], 
nanomagnets  [5-7],  mechanical  nanoresonators  [8],  and 
recently  in  quantum  measurements  [9-1 1]. 

The  analysis  of  switching  conventionally  relies  on  the 
assumption  that  the  underlying  noise  is  Gaussian.  Then  the 
switching  exponent  Q,  i.e.,  the  exponent  in  the  expression 
for  the  switching  rate  W  exp(—  Q),  displays  a  power-law 
dependence  on  the  distance  to  the  bifurcation  point  in  the 
parameter  space  77,  Q  77^  [1,12,13].  For  systems  in 
thermal  equilibrium,  /  =  3/2  for  a  saddle-node  bifurca¬ 
tion  and  /  =  2  for  a  pitchfork  bifurcation.  This  applies  also 
to  systems  far  from  equilibrium  [14,15]. 

Recently,  there  has  been  much  interest  in  large  fluctua¬ 
tions  and  switching  induced  by  non-Gaussian  noise  [lb- 
19],  Such  switching  can  be  used  to  determine  the  noise 
statistics  [20-26].  However,  the  features  of  the  switching 
rate  near  bifurcation  points  have  not  been  explored.  Yet, 
one  may  expect  that  the  77  dependence  of  the  switching 
exponent  will  differ  from  that  for  a  Gaussian  noise  and  will 
be  very  sensitive  to  the  noise  statistics. 

In  this  Letter  we  study  the  behavior  of  the  switching 
exponent  Q  for  systems  driven  by  Poisson  noise.  Such 
noise  is  often  encountered  in  photon  statistics  and  in  the 
statistics  of  current  in  mesoscopic  conductors.  We  show 
that  the  scaling  is  not  described  by  a  simple  power  law,  and 
the  overall  77  dependence  of  Q  is  much  weaker  than  for 


Gaussian  noise.  Surprisingly,  if  in  addition  to  Poisson  noise 
the  system  is  driven  even  by  a  comparatively  weak 
Gaussian  noise,  sufficiently  close  to  the  bifurcation  point 
this  noise  dominates  and  there  occurs  a  crossover  to  the 
standard  scaling  of  Q  for  Gaussian  noise. 

Generally,  one  would  expect  that,  unless  it  is  very  weak, 
a  Poisson  noise  would  make  a  stronger  effect  on  the 
switching  rate  than  a  Gaussian  noise.  Because  switching 
is  a  rare  event  on  the  scale  of  the  characteristic  relaxation 
time  of  the  system  tr,  it  requires  a  large  fluctuation,  whose 
probability  is  determined  by  the  tail  of  the  noise  distribu¬ 
tion.  Such  a  tail  is  less  steep  for  a  Poisson  noise  than  for  a 
Gaussian  noise. 

The  “takeover”  of  the  switching  rate  by  a  weak 
Gaussian  noise  close  to  a  bifurcation  point  is  a  more  subtle 
effect.  It  emerges  because  of  the  qualitatively  different 
ways  the  fluctuations  leading  to  switching  occur  for 
Gaussian  and  Poisson  noises.  This  can  be  understood 
from  the  equation  of  motion  for  the  slow  variable  q, 

q  =  -U'(q)  +  f(t).  (1) 

Here,  U{q)  is  the  effective  potential;  for  the  saddle-node 
and  pitchfork  bifurcations,  U  =  Uiss,(q)  and  U  = 
t/(PF>(<7),  respectively  [27],  with 

U(SN\q)  =  pq  -  ^<?3,  UiPF](q)  =  ~^VT  +  ^<?4- 

(2) 

For  77  >  0  the  system  has  a  stable  state  qa  at  the  minimum 
of  U(q)  (or  two  such  states,  for  a  pitchfork  bifurcation)  and 
a  saddle  point  qs  at  the  local  maximum  of  Uiq):  for  77  =  0 
these  states  merge  together.  The  results  can  be  immediately 
extended  also  to  the  case  [/(PF)  — > ►  —  f/(PF\  the  subcritical 
pitchfork  bifurcation  where  for  77  =  0  a  stable  state  merges 
with  two  unstable  states. 

The  force  f(t)  in  Eq.  (1)  is  noise.  We  will  consider  the 
cases  where  /(f)  is  a  Poisson  noise,  /(f)  =  fP{t)  with 
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fp{t)  =  gX„<5(f  —  tn),  or  a  white  Gaussian  noise,  f(t)  = 
fG(t),  =  2 DS(t  —  t'),  or  a  combination 

/(f)  =  fP{t)  +  /G ( t) .  Both  g  and  D  are  assumed  small, 
so  that  the  switching  rate  W  <C  t~l  =  U"(qa).  The  noises 
are  S  correlated  in  time,  because  the  system  motion  is 
slow;  they  are  also  independent  of  q,  since  of  interest  is  a 
small  region  in  the  system  phase  space  [14]. 

We  start  with  the  qualitative  picture  of  Poisson-noise- 
induced  switching.  Here  switching  occurs  only  for  the 
appropriate  pulse  polarity,  {qs  —  qa)/g  >  0.  A  single 
noise  pulse  shifts  the  coordinate  q  by  pulse  area  g. 
Switching  requires  driving  q  from  qa  to  qs,  from  where 
the  system  will  switch  with  probability  —1/2.  The  neces¬ 
sary  n0  —  (qs  —  qa)/g  pulses  should  occur  within  time 
£  tr,  so  that  the  system  cannot  relax  back  to  the  attractor 
between  the  pulses.  The  probability  of  such  pulse  sequence 
is  equal  to  (p/)"0  exp(—  vtr)/n0\,  where  v  is  the  average 
pulse  frequency.  By  construction,  this  is  the  probability  to 
switch  in  time  tr,  it  is  ~Wtr,  and  thus  gives  the  switching 
rate.  The  corresponding  estimate  of  the  switching  exponent 
Q  =  Q_P  for  n0  1,  vtr  is 

Q.p  ~Y(qs  ~  qa)/g\  H(<?s  ~  qa)/g vtr\  (3) 

From  Eq.  (2),  |g5  —  qa\  771/2,  whereas  tr  —  rjl~ £  with 
£  =  3/2  and  /  =  2  for  the  saddle-node  and  pitchfork 
bifurcations.  It  is  seen  from  Eq.  (3)  that  QP  77 1  2  but 
contains  a  large  77-dependent  logarithmic  factor;  the  ratio 
(qs  ~  qj/gvtr  is  the  large  parameter  of  the  theory. 

White  Gaussian  noise  fG(t),  on  the  other  hand,  leads  to 
switching  by  providing  a  force  that  overcomes  the  deter¬ 
ministic  force  U'(q).  From  Eq.  (2),  \U'(q)\  £  for 

qa  <  q  <  qs.  Since  the  probability  of  noise  realization  is 

exp[-  f  dtfc(t)/4D]  [28]  and  the  duration  of  the 
needed  noise  outburst  is  ~tr,  by  setting  this  probability 
to  be  ~Wtr  and  fG  —  U' ,  we  obtain  for  the  switching 
exponent  Q_  =  Q_G  the  familiar  expression  [13,14] 

Q.G  =  CV(/D  (4) 

(C  =  4/3  and  C  =  1/4  for  the  saddle-node  and  pitchfork 
bifurcation,  respectively,  see  below). 

With  decreasing  77,  the  distance  between  the  stationary 
states  qs  —  qa  77 1,72  decreases  slower  than  the  determi¬ 
nistic  force,  \U'(q)\  S  l/2.  Therefore  an  outburst  of 
Poisson  noise  required  for  a  transition  decreases  slower 
than  that  of  Gaussian  noise.  Respectively,  as  seen  from 
Eqs.  (3)  and  (4),  the  switching  exponent  decreases  much 
slower  for  Poisson  noise  than  for  Gaussian  noise.  As  a 
result,  even  for  small  Gaussian-noise  intensity  D  <C  g,  the 
switching  rate  is  determined  by  Gaussian  noise  for  suffi¬ 
ciently  small  17,  whereas  for  larger  97  it  is  determined  by 
Poisson  noise. 

The  crossover  between  Gaussian-  and  Poisson-noise- 
dominated  switching  is  clearly  seen  in  Fig.  1.  It  presents 
the  results  on  switching  near  a  saddle-node  bifurcation 
point  in  the  presence  of  both  noises.  Close  to  the  bifurca- 


In  t| 

FIG.  1  (color  online).  The  crossover  of  the  switching  exponent 
<2  from  the  Poisson-noise-  to  Gaussian-noise-dominated  behav¬ 
ior  with  the  decreasing  distance  77  to  the  saddle-node  bifurcation 
point.  The  solid  line  shows  the  result  of  Eqs.  (6)  and  (7)  for  the 
Gaussian-noise  intensity  D  =  0.009,  the  average  Poisson  pulse 
frequency  v  =  0.09,  and  the  Poisson  pulse  amplitude  g  =  0. 185; 
the  data  points  are  the  results  of  numerical  simulations  in  which 

was  found  from  W  using  the  Gaussian-noise  prefactor,  W  ~ 
(r)1 / tt)  exp(— Q) .  The  dashed  lines  and  the  dash-dotted  lines 
were  obtained  from  Eqs.  (6)  and  (7)  by  setting  D  =  0  or  g  =  0, 
respectively. 

tion  point  the  slope  of  lnl?  vs  In  77  is  3/2,  as  for  Gaussian 
noise.  However,  as  77  increases,  the  slope  approaches  that 
for  purely  Poisson  noise. 

We  now  give  a  quantitative  theory  of  the  switching 
exponent  and  find  its  77  dependence  for  Poisson  noise  as 
well  as  study  the  crossover  between  the  Poisson-  and 
Gaussian-noise  asymptotic.  As  a  starting  point,  we  use  a 
generalized  Fokker-Planck  equation  for  the  probability 
density  of  the  system,  which  follows  from  Eq.  (1): 

dtP  =  dq[U'(q)p(q)]  +  Dd2p(q)  +  v[p(q  -  g)  -  p(q)\ 

(5) 

The  last  two  terms  in  this  equation  describe  the  effect  of  the 
Poisson  noise,  i.e.,  of  uncorrelated  pulses  with  average 
frequency  v  that  shift  the  coordinate  by  g.  These  terms 
have  the  same  form  as  reaction  terms  in  the  master  equa¬ 
tion  for  a  reaction  system,  with  q  and  g  being  the  number 
of  species  and  the  change  of  this  number  in  a  reaction, 
respectively  [29],  except  that  in  the  present  case  q  is 
continuous;  the  analogy  applies  where  g  is  small  compared 
to  the  typical  scale  of  q ,  in  particular,  compared  to  qs  — 
qa.  In  this  case  a  Poisson  noise  in  the  Langevin  equation 
mimics  reactions  in  reaction  systems. 

For  small  D  and  g,  the  switching  rate  is  determined 
by  the  probability  current  away  from  the  initially  occupied 
attraction  basin  [30].  The  current  is  independent  of  time  for 
fr  «  t  «  W~ 1  and,  to  logarithmic  accuracy,  is  given  by 
the  quasistationary  probability  distribution  at  the  saddle 
point.  To  find  this  distribution,  we  solve  Eq.  (5)  in  the 
eikonal  approximation:  we  set  p(q)  =  exp[  — .v(<y)].  as¬ 
sume  that  s(q)  —  s(qa)  »  1,  and  keep  the  leading  order 
terms  in  dqs  [for  example,  we  disregard  d2s  compared  to 
( ft£/  v)2] -  This  leads  to  the  equation  for  .v ( <7 )  of  the  form 
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H{q,  dqs)  =  0,  where 

H(q,  p)  =  -pU'(q)  +  Dp2  -  v[\  -  exp(gp)].  (6) 

Equation  (6)  maps  the  problem  of  the  quasistationary 
distribution  of  the  fluctuating  system  onto  the  problem  of 
Hamiltonian  dynamics  of  an  auxiliary  conservative  system 
with  coordinate  q,  momentum  p,  and  Hamiltonian  H(q,  p). 
The  switching  exponent  is  Q_  =  s{qg)  —  s(qa),  or 

Q.  =  P  p(q)dq.  (7) 

J  9a 

Here,  qa,  qg  are  the  stationary  states  shifted  to  allow  for  the 
nonzero  mean  of  the  Poisson  noise,  cp  —  qi  ~  vg/U"{qj) 
with  i  =  a,  S.  The  momentum  p(q)  in  Eq.  (7)  is  the  non¬ 
trivial  solution  of  equation  H{q,  p)  =  0;  it  gives  the  tra¬ 
jectory  of  the  auxiliary  system  which  goes  from  qa  to  qs. 
Eqs.  (6)  and  (7)  could  be  obtained  also  by  finding  the 
probability  density  of  the  most  probable  realization  of 
noise  f(t)  necessary  to  drive  the  system  from  qa  to  q$ 
[19,23,31],  albeit  such  calculation  would  be  somewhat 
more  involved. 

For  purely  Poisson  noise,  i.e.,  where  D  =  0,  for  small 
pulse  area  g,  the  equation  H{q,  p)  =  0  gives  p  ~ 
g_1{lnu(g)  +  ln[lnu(g)]},  where  v\q)  =  U'{q)/ vg.  This 
estimate  applies  provided  v{q) »  1;  we  disregarded 
higher-order  corrections  “T/ln [u(g)].  The  condition 
v{q)  »  1  holds  in  much  of  the  region  between  qa  and 
q§  except  the  immediate  vicinities  of  qa,  qs,  i.e.,  for 
I q  -  qa\,  \q  -  qs\  »  v\g\/U"{qa).  This  follows  from 
the  estimate  \U'{q)\  ;£  17 1/2  in  the  central  part  of  the 
interval  ( qa ,  qs)  and  the  inequality  {qs  ~  qa)/§  ^  vtr  = 
v/U"(qa)  discussed  above. 

Keeping  in  p{q)  the  leading  order  term  and  replacing 
U'{q)  by  its  maximal  value  between  qa  and  qg,  from 
Eq.  (7)  we  obtain  the  following  estimates  for  the  switching 
exponent  for  the  saddle-node  and  pitchfork  bifurcation, 
respectively, 

QpN)  ~  {2px/2  /  g)\n{Kim  p  /  vg), 

£pPF)  ~  (pl^2/g)  ln(«-(PF)p3/2/pg). 

The  parameters  /c(SN),  Ktw)  in  the  arguments  of  the  loga¬ 
rithms  are  ~1.  A  simple  choice  /<(SN)  =  2  and  /c(PF)  =  1 
gives  a  close  agreement  of  Eq.  (8)  with  the  results  obtained 
by  numerically  solving  equation  H  =  0  and  Eq.  (7),  which 
are  shown  in  Figs.  2  and  3.  For  the  chosen  v,  the  difference 
is  <10%  for  £  20  and  g  >0.1.  Equation  (8)  justifies 
the  estimate  Eq.  (3). 

The  power-law  factor  in  Eq.  (8)  is  the  same  for  both 
types  of  the  bifurcation  points  [32].  It  is  determined  simply 
by  the  distance  between  qs  and  qa.  However,  the  argu¬ 
ments  of  the  logarithms  are  different.  The  logarithmic 
factors  significantly  change  the  switching  exponent  com- 


FIG.  2  (color  online).  The  switching  exponent  for  a 

Poisson  noise  as  a  function  of  the  distance  17  to  the  bifurcation 
point  for  the  saddle-node  bifurcation.  The  mean  frequency  of 
noise  pulses  is  v  =  0. 1 .  The  solid  lines  are  obtained  from  Eq.  (7) 
using  a  numerical  solution  of  equation  H(q,  p)  =  0.  The  data 
points  show  the  results  of  numerical  simulations  of  switching; 
the  plotted  quantity  is  lnfp'/z /Wtt). 

pared  to  a  simple  power-law  scaling  QP  p]'2 / g-  For 
example,  for  the  pitchfork  bifurcation  for  g  =  0.2  in  Fig.  3 
the  logarithmic  factor  varies  from  3.6  for  p  =  0.5  to  2.2  for 

p  =  0.2. 

The  results  of  the  asymptotic  theory  were  compared 
with  numerical  simulations  of  the  system  Eq.  (1).  As 
seen  from  Figs.  2  and  3,  they  are  in  excellent  agreement. 
In  determining  the  switching  exponent  from  the  switching 
rate  W,  we  used  the  prefactor  in  W  which  coincides  with 
that  for  white  Gaussian  noise  [33].  We  note  that,  for  the 
supercritical  pitchfork  bifurcation  described  by  the  poten¬ 
tial  (/<PF)(g)  in  Eq.  (1),  a  unipolar  (all  pulses  of  the  same 
sign)  Poisson  noise  leads  to  switching  from  only  one  of  the 
two  coexisting  stable  states  at  the  minima  of  (/(PF)(<?)-  For 
the  subcritical  bifurcation  described  by  the  potential 
—  l!>PF)(q),  Poisson  noise  always  leads  to  decay  of  a  meta¬ 
stable  state. 


FIG.  3  (color  online).  The  switching  exponent  for  a 

Poisson  noise  as  a  function  of  the  distance  p  to  the  bifurcation 
point  for  the  pitchfork  bifurcation.  The  mean  frequency  of  noise 
pulses  is  v  =  0.05.  The  solid  lines  are  obtained  from  a  numerical 
solution  of  equation  H(q,  p)  =  0  followed  by  integration,  Eq.  (7). 
The  data  points  show  the  results  of  numerical  simulations  of 
switching;  the  plotted  quantity  is  ln( p/Wir^). 
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In  the  opposite  case  of  purely  Gaussian  noise,  g  =  0,  the 
nontrivial  solution  of  equation  H{q,  p)  =  0  is  p  = 
U'(q ) /£).  It  leads  to  a  power-law  dependence  of  the  switch¬ 
ing  exponent  on  p  described  by  Eq.  (4)  and  to  the  appro¬ 
priate  values  of  the  constant  C  in  this  equation. 

Equation  (6)  describes  quantitatively  the  crossover  from 
Poisson-  to  Gaussian-noise-dominated  switching  as  p  ap¬ 
proaches  the  bifurcation  value  77  =  0.  If  the  Gaussian  noise 
is  much  weaker  than  the  Poisson  noise,  D  <C  g,  far  from 
qa,  qs  the  momentum  p{q)  is  close  to  the  Poisson  noise 
(D  =  0)  solution  provided  p  is  not  too  small,  so  that 
M{p)  »  D/g,  with  M(p)  =  p£~ll2 /  \n(p€~1/2 / vg)  (we 
use  that,  near  its  maximum,  \  JJ'\  ~  p^  l/2).  On  the  other 
hand,  where  M(p)  <5C  D/g  the  momentum  is  determined 
by  the  Gaussian-noise  (g  =  0)  solution.  The  position  of  the 
crossover  on  the  p  axis  is  given  by  M{p)  ~  D/g.  Not 
surprisingly,  for  such  p  the  switching  exponents  for  the 
purely  Poisson  and  purely  Gaussian  noises  become  of  the 
same  order  of  magnitude,  Q_P  ~  Q_G.  This  argument  is 
confirmed  by  the  data  in  Fig.  1. 

The  situation  where  a  weak  Gaussian  noise  is  present 
even  where  other  sources  of  noise  are  dominating  is  typical 
for  practically  any  physical  system.  Such  noise  very  often 
comes  simply  from  the  coupling  of  the  system  to  a  reser¬ 
voir  that  leads  to  energy  dissipation.  The  results  of  this 
Letter  explain  why  near  bifurcation  points  there  is  often 
observed  the  power-law  scaling  typical  for  Gaussian  noise 
even  where  this  noise  is  comparatively  weak. 

In  conclusion,  we  have  considered  noise-induced 
switching  due  to  a  non-Gaussian  noise  near  two  generic 
types  of  the  bifurcation  points:  saddle  node  and  pitchfork. 
In  contrast  to  the  case  of  Gaussian  noise,  where  the  switch¬ 
ing  exponent  scales  as  a  power  of  the  distance  to  the 
bifurcation  point  17,  for  a  non-Gaussian  noise  the  exponent 
generally  displays  a  more  complicated  dependence  on  p. 
We  have  found  it  for  a  Poisson  noise,  in  which  case,  along 
with  a  power-law  factor,  the  exponent  has  a  large  logarith¬ 
mic  factor.  It  turned  out  that  even  a  weak  additional 
Gaussian  noise  becomes  the  major  cause  of  switching 
sufficiently  close  to  the  bifurcation  point.  A  qualitative 
and  quantitative  description  of  the  crossover  from 
Poisson-  to  Gaussian-noise-controlled  switching  and  of 
the  p  dependence  of  the  switching  exponent  is  in  full 
agreement  with  numerical  simulations. 
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